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Abstract. We present detailed theoretical mass-radius relations for massive white dwarf stars with oxygen-neon 
cores. This work is motivated by recent observational evidence about the existence of white dwarf stars with very 
high surface gravities. Our results are based on evolutionary calculations that take into account the chemical 
composition expected from the evolutionary history of massive white dwarf progenitors. We present theoretical 
mass-radius relations for stellar mass values ranging from 1.06 to 1.30 Mq with a step of 0.02 Mq and effective 
temperatures from 150 000 K to « 5 000 K. A novel aspect predicted by our calculations is that the mass-radius 
relation for the most massive white dwarfs exhibits a marked dependence on the neutrino luminosity. Extensive 
tabulations for massive white dwarfs, accessible from our web site, are presented as well. 
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1. Introduction 

White dwarfs are the end-product of the evolution of 
low- and intermediate-mass stars. Thus, they preserve im- 
portant clues about the formation and evolution of our 
Galaxy. This information can be retrieved by studying 
their observed mass, kinematic and luminosity distribu- 
tions, provided that we have good structural and evolu- 
tionary models for the progenitors of white dwarfs and 
for the white dwarfs themselves. In particular, two funda- 
mental tools to analyze the properties of the white dwarf 
population as a whole are the white dwarf luminosity func- 
tion and the mass distribution of white dwarfs. The former 
has been consistently used to obtain estimates of the age of 
the Galactic disk (Winget et al. 1987; Garcfa-Berro et al. 
1988; Hernanz et al. 1994; Richer et al. 2000) and the past 
history of the galactic star formation rate (Noh & Scalo 
1990; Dfaz-Pinto et al. 1994). The latter yields important 
information about the late stages of stellar evolution since 
it reveals the amount of mass lost during stellar evolution 
from an initial stellar mass distribution (Licbert, Bergeron 
& Holberg 2005). Moreover, white dwarfs have been pro- 
posed to contribute significantly to the mass budget of 
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our Galaxy — see Garci'a-Berro et al. (2004) for a recent 
discussion on this topic. 

However, the fundamental parameters of white dwarf 
stars, such as their individual masses and radii, are dif- 
ficult to estimate. Hence, reliable mass and radius deter- 
minations are only available for a relatively small fraction 
of the whole local white dwarf population. Consequently, 
one of the most important theoretical results, namely the 
mass-radius relation, is poorly constrained by the observa- 
tions, despite the improved evolutionary sequences avail- 
able nowadays (Benvenuto & Althaus 1999; Hansen 1999; 
Salaris et al. 2000). 

A good determination of the white dwarf mass dis- 
tribution provides useful constraints to the initial-final 
mass relationship, a key ingredient for modelling the 
different white dwarf populations (thin and thick disk, 
spheroid...), and a touchstone for the stellar evolution- 
ary models. Several works have focused on this issue. 
For instance, Finley et al. (1997) determined the effec- 
tive temperatures and surface gravities of 174 DA white 
dwarfs with Tcff > 25 000 K, and from that they obtained 
the mass distribution of hot DA white dwarfs. Marsh et 
al. (1997) describe the results for a sample of 89 EUV- 
selected white dwarfs and Vennes (1999) determined the 
mass distribution for a sample of 141 EUV-selected white 
dwarfs. Silvestri et al. (2001) studied the gravitational red- 
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shifts and mass distribution of 41 white dwarfs in com- 
mon proper motion binary systems. Finally, Liebert et al. 
(2005) studied a complete sample of 348 DA white dwarfs 
from the Palomar Green survey. All these studies report 
the existence of a clear, well defined peak at ~ 0.6 Mq 
with a typical dispersion of ^ 0.1 Mq. Additionally, a sec- 
ondary peak is also reported at masses between 1.0 and 
1.2 Mq. These high-mass objects should have ONe cores, 
the primary products of the ^^C nuclear reactions, in con- 
trast with average-mass white dwarfs, which undoubtedly 
have CO cores. The origin of these high-mass white dwarfs 
is uncertain, and it has been suggested that they could be 
the result of either binary evolution (Marsh et al. 1997) 
— but see Guerrero et al. (2004) — or of the evolution of 
heavy-weight intermediate-mass single stars (Ritossa ct al. 
1996; Garcia-Berro et al. 1997; Iben et al. 1997, Ritossa 
et al. 1999). In any case these white dwarfs do exist and 
a good mass-radius relation is needed in order to accu- 
rately determine their respective masses. All the mass de- 
terminations, including the most recent ones, rely on the 
evolutionary sequences of Wood (1995), which are valid 
for CO cores. To the best of our knowledge no reliable 
mass-radius relations for massive white dwarfs with ONe 
cores exist today and, thus, this is precisely the goal of the 
present work. Our paper is organized as follows. In §2 we 
present our input physics. In §3 we discuss our evolution- 
ary sequences. Finally, in the last section we summarize 
our findings and draw our conclusions. 

2. Input physics 

The evolutionary code used in this paper is that used 
in our recent work about the evolution and pulsational 
properties of massive white dwarfs with oxygen and neon 
cores (Corsico et al. 2004), and we refer the reader to 
that paper for details about the input physics and about 
the procedure followed to obtain the initial white dwarf 
configurations. The code is based on updated and de- 
lailed constitutive physics. Briefly, the equation of state 
for the low-density regime includes an updated version of 
the equation of state of Magni & Mazzitelli (1979). For 
the high-density regime, partially degenerate electrons, 
radiation pressure, ionic contributions and Coulomb in- 
teractions arc considered. Radiative opacities are those 
of OPAL (Iglesias & Rogers 1996) including carbon- and 
oxygen-rich compositions, complemented at low temper- 
atures with the Alexander & Ferguson (1994) molecular 
opacities. Neutrino emission rates for pair, photo, plasma 
and bremsstrahlung processes, and high-density conduc- 
tive opacities are taken from Itoh et al. (1996a,b), and Itoh 
et al. (1994). Convection is treated in the framework of the 
mixing length theory as given by the ML2 parameteriza- 
tion (Tassoul et al. 1990). In this paper, we have neglected 
the effects of crystallization. Because we are interested in 
the calculation of mass radius relations, this assumption 
bears virtually no relevance. However, it is worth noting 
that crystallization in massive white dwarfs takes place 
at such high stellar luminosities that its impact on the 
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Fig. 1. Hertzsprung- Russell diagram for two 1.28 Mq evo- 
lutionary sequences with different artificial heating. Lines 
of constant radius for R = 0.01, 0.001 and 0.0001% 
(from top to bottom) are also shown. See text for details. 

white dwarf cooling is less relevant. In addition, chemical 
redistribution due to phase separation upon crystalliza- 
ton of the oxygen-neon core is of minor importance — 
see Corsico et al. (2004) for details — and has not been 
considered in this work. 

Our starting stellar configurations correspond to hot 
white dwarf structures with a realistic chemical stratifi- 
cation appropriate to white dwarfs resulting from progen- 
itor stars of solar metallicity that are expected to have 
burnt carbon in scmidegenerate conditions. All the evolu- 
tionary sequences considered in this work have the same 
chemical abundance profiles, and correspond to that illus- 
trated in figure 4 of Corsico et al. (2004). Although minor 
changes in the chemical profile are expected because of 
the different masses of progenitor objects, these have a 
negligible influence in the mass-radius relation. The outer 
layer chemical stratification consists of a pure hydrogen 
envelope of 1.4 x 10"®M, overlying a pure helimn shell of 
4 X lO-^M* and, below that, a buffer rich in and ^^'O. 
We mention that the amount of hydrogen we adopted is 
an upper limit as imposed by nuclear reactions. We also 
find that the mass-radius relations are almost insensitive 
to the hydrogen envelope mass. The core is primarily com- 
posed of ^^O and ^°Ne, plus some traces of ^^C, ^^Na and 
^''Mg. The chemical composition of the core is the result 
of repeated shell flashes that take place during the carbon 
burning phase in massive white dwarf progenitors (Garcia- 
Berro et al. 1997). The oxygen-neon core of our models is 
characterized by flat profiles that arc the result of a reho- 
mogenization process induced by Ray leigh- Taylor instabil- 
ities in regions with negative molecular weight gradients. 
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Fig. 2. Internal temperature in terms of the outer mass 
fraction. The top and bottom panel correspond, respec- 
tively, to three selected models from the evolutionary se- 
quences depicted with thin and thick lines in figure 1. 
Solid, dashed and dot-dashed lines correspond to lumi- 
nosities of about logL/LQ =0.49, —0.16 and —0.58, re- 
spectively. Note that for the two most luminous models, 
the internal temperature is strongly influenced by neutrino 
emission. 

a process that is expected to take place along the white 
dwarf evolutionary sequence. We mention that the lowest 
stellar mass we considered in this study is close to the 
theoretical lower limit of w 1.05 A/q for an oxygen-neon 
white dwarf to be formed (Gil-Pons et al. 2003). 

The initial stellar models needed to start our cool- 
ing sequences were obtained by means of the same ar- 
tificial evolutionary procedure described in Gautschy & 
Ahhaus (2002) and Ahhaus & Corsico (2004). Briefly, we 
start from a hot 0.95 Mq white dwarf model obtained self- 
consistently by evolving the progenitor star initially on the 
ZAMS (Althaus et al. 2003) to the white dwarf regime. In 
order to obtain initial white dwarf configurations for the 
different masses studied here we scaled the stellar mass 
and incorporated an artificial energy generation rate un- 
til the white dwarf reaches a luminosity far larger than 
that of the initial models considered here as meaning- 
ful. We checked the thermal structure of the models and 
found that the transitory state has negligible influence on 
the thermal conflguration of the white dwarf interior well 
before the initial configuration is reached (see next sec- 
tion) . Specifically, we checked the thermal structure of the 
1.28 Mq model by extending the heating artificial proce- 
dure even further beyond the knee at high luminosities in 
the Hertzsprung- Russell diagram, as shown in Fig ^ In 
this figure, the thick line displays the evolutionary track 
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Fig. 3. Internal temperature in terms of the outer mass 
fraction. Solid, dashed and dot-dashed lines correspond 
to the 1.06 Mq white dwarf models at logL/LQ =0.66, 
—0.02 and —0.38, respectively. Note that for the most lu- 
minous model, the maximum temperature is off-center. 

for one of our standard sequences (that of the 1.28 Mq 
model), whilst the thin line corresponds to the cooling 
track obtained when the artificial heating is extended to 
much larger luminosities. As can be seen the two evolu- 
tionary tracks are indistinguishable, thus confirming that 
the evolutionary results presented below do not depend 
on the initial conditions. 

We have computed the evolution of massive white 
dwarf configurations from hot effective temperature stages 
down to very low luminosities. The values of the stellar 
mass considered range from 1.06 to 1.30 Mq with a step 
of 0.02 Mq. The lower limit of this mass range corresponds 
to the minimum mass for a white dwarf to have an ONe 
core, whereas the upper limit corresponds to white dwarf 
masses for which post-newtonian corrections start to be 
important. Nevertheless, it is worth noting that extremely 
massive white dwarfs — those with M* ^ 1.3 Mq — are 
very unusual. The mass-radius relations are presented for 
effective temperature values ranging from about 150 000 K 
to « 5 000 K. 

3. Evolutionary results 

We begin by examining Fig. |21 which shows the inter- 
nal temperature as a function of the outer mass fraction, 
log(l — Mr/M*), for three selected evolutionary stages of 
the 1.28 A/q white dwarf sequence corresponding, from 
top to bottom, to total luminosities of about logi/^Q = 
0.49, —0.16 and —0.58, respectively, which are marked as 
dots in Fig. 1 . As can be seen the three models are already 
at advanced evolutionary stages, and, consequently, their 
thermal structure should be reliable. The top panel refers 
to our standard sequence, whilst the bottom panel dis- 
plays the situation for models belonging to the sequence 
that started at very high luminosities (denoted by thin line 
in Fig. 1). Note that the temperature profiles of the models 
are the same for both sequences, reflecting the correctness 
of our artificial procedure to generate starting white dwarf 
configurations. 



4 



Althaus et al.: Mass-radius relations for massive white dwarf stars 




o -i 




(J 

o -i 



M= 1.06 M 




I I I I L^r I I I I I I 

0.2 0.4 0.6 0.8 1 

MM 

r 

Fig. 4. Neutrino emission rate (in erg/g/s) in terms of 
the mass fraction for three selected 1.28 Mq stellar mod- 
els. Solid, dotted and dashed lines correspond to the to- 
tal rate of neutrino emission, and to the rates due to 
brehmsstrahlung and plasma processes, respectively. The 
top, central and bottom panel correspond to luminosities 
of about log L/Lq =0.51, -0.17 and -0.58. 

It is well known that neutrino losses considerably af- 
fect both the cooling time-scales and the structural prop- 
erties of hot white dwarfs. Indeed, neutrino emission con- 
stitutes the most important physical process affecting the 
evolution of very hot CO white dwarfs, regardless of their 
mass. As can be seen in Figs. |21 and 13 this statement is 
also true for the massive ONe white dwarfs studied here. 
Note that, for the most massive model and for the lumi- 
nosity range considered here, neutrino losses do not cause 
the maximum temperature to occur off-center (top panel 
of Fig. El , as it is the case for the less massive model of 
1.06 Mq, which is displayed in Fig. O In particular it is 
worth noting that for the 1.28 Mq white dwarf cooling se- 
quence, the model with the largest photon luminosity — 
solid lines in Fig. El — temperature shows a local maxi- 
mum at log(l — Mi/Af*) ~ —2.2, well inside the ONe core. 
This can be understood by examining Fig.01 which shows 
the neutrino emission rates (in erg/g/s) as a function of 
the mass fraction for the 1.28 A/q white dwarf model se- 
quence. Note that the maximum neutrino emissivity oc- 
curs towards the outer edge of the degenerate core, with 
plasma neutrino being the dominant contribution. Indeed, 
for massive white dwarfs, the emission rate resulting from 
plasma neutrino is characterized by a maximum at den- 
sities much lower than the central densities of such mod- 
els, as clearly shown, for instance, in Fig. 6 of Itoh ct al. 
(1996b). 



-1.5 

Fig. 5. Neutrino emission luminosity resulting from 
plasma and bremsstrahlung processes (solid and dashed 
lines, respectively) for the 1.06 and 1.28 Mq white dwarf 
sequences. In the interests of comparison, we depict the 
situation for which neutrino energy losses are equal to the 
photon luminosity of the white dwarf. 



During the first evolutionary stages of the white dwarf 
models studied here — those with log L/Lq ^ —0.75 — 
most of the neutrino emission proceeds via the plasma 
process, as it is apparent from Fig. |31 (see also Fig. 
However, for the most massive white dwarfs, neutrino 
losses through the bremsstrahlung process are not negligi- 
ble whatsoever, being comparable to or even larger than 
the photon luminosity of the white dwarf. This is due to 
the strong dependence of the neutrino emissivity on the 
density and on the average charge of the ions. Moreover, 
the neutrino luminosity far exceeds the photon luminosity 
during these early stages of the evolution. 

Neutrino losses constitute a powerful sink of energy 
in hot massive white dwarfs and leave their signatures 
in the mass-radius relations of the sequences, as it is ap- 
parent from Fig. |B| Indeed, as neutrino emission gradu- 
ally decreases the models undergo an appreciable change 
in their radius. This effect is clearly noticeable in the 
1.28 A/q white dwarf model shown in the bottom panel 
of Fig.|Sl which shows a marked reduction in the radius of 
the model as the neutrino emission becomes extinct below 
logL/LQ w —0.25. Note that this reduction in the stellar 
radius is absent if neutrino losses are not considered (the 
corresponding evolutionary sequence is represented by a 
dashed line in the botton panel of Fig.|SJ). This coupling, 
which is not found for stellar masses smaller than about 
1.20 A'/q, has not been reported before. We stress that 
this behaviour is only noticeable for very massive white 
dwarf models, because in these models the maximum of 
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Fig. 6. Upper panel: Evolution of the neutrino emis- 
sion luminosity in terms of the stellar luminosity for the 
1.28Mq white dwarf sequence. Bottom panel: stellar ra- 
dius as a function of stellar luminosity for the 1.28 Mq 
white dwarf sequence. The dashed line corresponds to the 
case in which neutrino losses have not been considered in 
the evolutionary sequence. Note the change in the stellar 
radius at log L/Lq « —0.25 as neutrino losses fade away. 



neutrino emission occurs in the outer, less degenerate re- 
gions of the core. On the contrary, in less massive white 
dwarfs, neutrino energy losses are large in the degener- 
ate core of the white dwarf. The upper panel of Fig. 
illustrates the behaviour of neutrino luminosity, L^, as a 
function of photon luminosity. Here, the various neutrino 
processes are shown in the same plot. Note that neutrino 
energy losses in massive white dwarfs far exceed the pho- 
ton luminosity. As the evolution proceeds, the neutrino 
luminosity fades away faster than the photon luminosity 
and neutrino cooling has little effect on the subsequent 
evolution of the models. This behaviour can be explained 
with the help of Fig. [Sjas follows. As it has been already 
shown, for the most massive ONe white dwarfs the neu- 
trino emissivity peaks in the outer partially degenerate 
layers of the white dwarf. Hence, these layers are cooler 
than in the case in which neutrino emission has been dis- 
regarded. Consequently, the radius of the white dwarf is 
smaller for the sequence in which neutrino emission has 
been properly taken into account. Since the peak of the 
neutrino emissivity is located at the edge of the degener- 
ate core, in a region in which the temperature gradient 
is large, the nearly isothermal core provides enough en- 
ergy to prevent the temperature to drop very rapidly in 
these regions. As a consequence the rate of change of the 
stellar radius is smaller than in the case in which neu- 
trinos are disregarded. Moreover, since these layers are 



' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' 






5.2 5.0 4.8 4.6 4.4 4.2 4.0 3.8 



Fig. 7. Surface gravity as a function of effective tempera- 
ture for all of the massive white dwarf evolutionary se- 
quences considered in this work. From bottom to top, 
curves correspond to sequences with stellar masses from 
1.06 to 1.30 Mq with a mass step of 0.02 Mq. As neu- 
trino emission becomes extinct at high effective tempera- 
tures, models becomes more compact, as it is reflected by 
the sudden change in the surface gravity, particularly for 
most massive models. For such sequences and advanced 
stages of evolution, finite-temperature effects are barely 
noticeable. 



already cooler than those of the model in which neutrino 
emission is disregarded the rate of change of the radius is 
smaller as well. Nevertheless, in these layers neutrino cool- 
ing is efficient and, therefore, the temperature decreases 
and the density increases. This, ultimately leads to a de- 
crease of the neutrino luminosity, which finally drops be- 
low the photon luminosity. At this stage the radius of the 
evolutionary sequence in which neutrinos are taken into 
account is indistinguishable from that in which neutrinos 
has been disregarded. 

The main result of our evolutionary calculations is 
summarized in Fig. 13 which shows the evolution of the 
surface gravity (in cgs units) for all our sequences as a 
function of the effective temperature. As it is well known, 
for a given effective temperature, more massive models 
are characterized by smaller radii and larger gravities. As 
cooling gradually proceeds, the surface gravity becomes 
larger, eventually reaching an almost constant value cor- 
responding to that of the zero-temperature configuration, 
as it should be expected for strongly degenerate configura- 
tions. Note that finite-temperature effects are appreciable 
only at high effective temperatures, in the early evolution- 
ary phases. 
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Table 1. Surface gravities for white dwarf models with 
different core chemical composition at T^g — 20 000K. 
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Fig. 8. Surface gravity as a function of effective tempera- 
ture for the 1.30, 1.24 and 1.10 Mq white dwarf sequences 
with various core chemical compositions. Solid, dashed 
and dot-dashed lines correspond to core compositions of 
oxygen-neon, oxygen and carbon, respectively. Note the 
effect of the core chemical composition on surface gravity, 
particulary for very massive white dwarfs. 



Because of the presence of abundant ^"Ne in the core 
of our massive white dwarf models, gravities are larger 
than those obtained when a CO core is adopted. This is 
a result of the non-ideal, corrective terms of the equa- 
tion of state. In fact, the presence of ^°Ne yields a much 
stronger interacting plasma than that of a regular CO mix- 
ture adopted in the previous works (Wood 1995). This 
is borne out by Fig. |S| which illustrates the evolution of 
the surface gravity for selected massive white dwarf mod- 
els with different core chemical composition. Specifically, 
solid, dashed and dot-dashed lines correspond to core com- 
positions of oxygen-neon, oxygen and carbon, respectively. 
As expected, for a given stellar mass, surface gravities are 
larger the higher the atomic number Z of the chemical con- 
stituent. Interestingly enough, very massive white dwarfs 
are rather sensitive to the adopted core chemical composi- 
tion. Such white dwarfs are characterized by densities high 
enough for electrons to be extremely relativistic through- 
out most of their interiors. Hence, Coulomb interactions 
grow in importance as the stellar mass increases. Finally, 
in table 1 we list the surface gravity values expected at 
Toff = 20 000K for some selected white dwarf sequences 
with different core chemical composition. 

4. Conclusions 

We have computed reliable mass-radius relations — or, 
equivalently, logg/Toff relations — for massive white 
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dwarfs with degenerate ONe cores. Our calculations en- 
compass masses ranging from 1.06 to 1.30 Mq, the ex- 
pected range of masses for which ONe white dwarfs should 
presumably exist. To the best of our knowledge this is 
the first attemtpt to compute such a relation with a re- 
alistic equation of state — which includes the non-ideal, 
corrective terms and the full temperature dependence — 
and reliable chemical profiles for the degenerate interior 
reflecting the previous evolutionary history of the pro- 
genitors of such white dwarfs. The results presented here 
will undoubtedly be helpful in the interpretation of recent 
observations of white dwarf stars with very high surface 
gravities (Dahn et al. 2004; Madej et al. 2004; Nalezyty & 
Madej 2004), which up to now relies on prior evolution- 
ary sequences which were computed assuming CO cores. 
Moreover, we have shown that the mass-radius relation 
shows a marked dependence on the neutrino luminosity, 
in such a way that for the most massive stars within this 
mass range a sudden decrease in the mass-radius relation- 
ship shows up shortly before the neutrino luminosity drops 
below the photon luminosity. 

Finally, we have prepared detailed tabulations of 
radius and gravity for all our massive white dwarf 
sequences, which are available at our web site: 
http : / / www . f caglp . unlp . edu . ar/evolgroup/. 
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